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Here, we propose a scheme to generate a controllable Ising interaction between superconducting 
flux qubits. Existing schemes rely on inducting couplings to realize Ising interactions between flux 
qubits, and the interaction strength is controlled by an applied magnetic field On the other hand, we 
have found a way to generate an interaction between the flux qubits via capacitive couplings. This 
has an advantage in individual addressability, because we can control the interaction strength by 
changing an applied voltage that can be easily localized. This is a crucial step toward the realizing 
superconducting flux qubit quantum computation. 


I. INTRODUCTION 

To realize fault-tolerant quantum computation, it is 
crucial to investigate a scheme to generate a cluster state 
in a scalable way. The cluster state is a universal resource 
for quantum computation, and this state can be used 
for a fault-tolerant scheme such as a surface code and 
topological code. One can generate a cluster state if we 
can turn on/off an Ising type interaction between qubits. 

Superconducting circuit is one of the promising sys¬ 
tems to realize such a cluster-state quantum computa¬ 
tion. Josephson junctions in the superconducting circuit 
can induce a non-linearity, and so one can construct a 
two-level system. There are several types of Josephson 
junction qubit: charge qubit 0, superconducting spin 
qubit superconducting flux qubit 0-0 , supercon¬ 
ducting phase qubit ISl [l(l|. superconducting transmon 
qubit fill . HI , fluxonium qubit 0111, and several hy¬ 
brid systems 0 , 01 ■ 

The transmon qubit 000, which is a cooper-pair 
box and relatively insensitive to low-frequency charge 
noise, is considered one of the powerful method of the 
qubit implementation by using superconducting circuit. 
Scheme of the tunable qubit- qub it capacitive coupling is 
proposed and demonstrated 0-0. The high fidelity 
qubit readout using a microwave amplifier is demon¬ 
strated 0-0]- Furthermore, high fidelity (99.4%) two- 
qubit gate using five qubits system is achieved. This 
result is the first step toward surface code scheme 0 . 
These results show a good scalability towards the real¬ 
ization of generating a large scale cluster state. 

The flux qubit consist of a superconducting loop con¬ 
taining several Josephson junctions. This system has a 
large anharmonicity and can be well approximated to a 
two-level system. Single qubit gate operations can be re¬ 
alized with high speed and reasonable fidelity 0. Mean- 
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while, the best observed coherence time is an order of 10 
pts 0,0. Furthermore, the tunable coupling schemes 
for two qu bit g ate operations are proposed and demon¬ 
strated |27 h 35| . Quantum non-demolition measurement 
of flux qubit during the coherence time is realized by us¬ 
ing Josephson bifurcation amplifier 0 - 0 ]. 

There are two typical tunable qubit-qubit coupling 
schemes, inductive coupling and capacitive coupling. In 
flux qubit system, existing schemes rely on inductive cou¬ 
pling with the external magnetic field. Several schemes of 
the tunable qubit- qub it inductive coupling are proposed 
and demonstrated 0 - 0 . However, it is hard to apply 
magnetic field to a localized region. Due to this prop¬ 
erty, it is difficult to achieve individual addressability of 
all qubits, because magnetic field may affect not only the 
target qubits but also other qubits as well. Therefore, 
it is important to perform two-qubit gates without af¬ 
fecting other qubits by using localized fields for scalable 
quantum computation. 

Here, we propose a way to generate and control the 
Ising type interaction between four-junction flux qubits 
using capacitive coupling. By using an applied voltage, 
we control the interaction between flux qubits that are 
connected by capacitance. Unlike the standard schemes, 
our scheme does not require to change the applied mag¬ 
netic field on the flux qubit for the control of the interac¬ 
tion. This may have advantage to implement two-qubit 
gates on the target qubits without affecting other qubits 
because applying local voltages is typically much easier 
than applying local magnetic flux. We take into account 
of realistic noise on this type of flux qubits, and estimate 
a qubit-parameter range where one can perform fault- 
tolerant quantum computation. Furthermore, we show a 
way to generate a two dimensional cluster state in a scal¬ 
able way. Our result paves the way to achieve the scalable 
quantum computation with superconducting flux qubits. 

The rest of this paper is organized as follows: In Sec¬ 
tion , we presents the design details of our flux qubit and 
effects on a flux qubit from the change of the parameters. 
In Section Hill we propose our scheme for generating Ising 
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type interaction between capacitively coupled supercon¬ 
ducting flux qubits. Moreover, we show the relationship 
between coupling strength and two types of errors caused 
by operation accuracy, the fluctuation of applied voltage 
and timing jitter. In Section m we present the analysis 
of our scheme for use in multi-qubit system. Addition¬ 
ally, we discuss how to suppress the non-nearest neighbor 
interactions by changing parameters and performing tt 
pulses. Furthermore, we show our procedure for generat¬ 
ing a one and two-dimensional cluster state using qubits 
on square lattice in less time. 


II. VOLTAGE CONTROLLED a-TUNABLE 
FLUX QUBIT 


Let us first show the circuit of a flux qubit that we pro¬ 
pose in Fig. 1(a) Here, X-shaped crosses denote Joseph- 



(a) (b) 


FIG. 1. (a)The circuit of a flux qubit in our design. This flux 
qubit has four Josephson junctions (JJ). Ej(n) and Cj( n ) de¬ 
note the Josephson energy and capacitance of n th Josephson 
junction JJn. The loop is threaded by an external magnetic 
flux /, and we can control the energy bias of the qubit via 
the magnetic flux. Node 1 represents the superconducting 
island. The electric potential of node 1 is Vi. (b) The a de¬ 
pendence of Eoi and E 12 where £01 denotes energy difference 
between the ground state and the first excited state, £12 de¬ 
notes energy difference between the first excited state and the 
second excited state. Here, we set £,( 1 ) = Ej ( 4 ) = 200 GHz, 
Ej( 2 ) = Ej ( 3 ) = 40 GHz, and £j(fc)/£ c (fc) = 80(A; = 1,2, ..,4). 


son junctions (JJ). The first Josephson junctions (JJ1) 
and the fourth Josephson junctions (JJ4) both have the 
same Josephson energies Ej and capacitances Cj. The 
second Josephson junctions (JJ2) and the third Joseph¬ 
son junction (JJ3) both have the same Josephson ener¬ 
gies and capacitances that are a times larger than those 
of JJ1 and JJ4. Josephson phases <p n , which is given by 
the gauge-invariant phase of each JJn, are subject to the 
following equation: 


(1) 


due to fluxoid quantization around the loop containing 
phases of Josephson junctions. / denotes the external 
magnetic flux through the loop of the qubit in units of 
the magnetic flux quantum $ 0 = ^ ■ The total Josephson 
energy U can be described as follows: 
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u = E j(k) (! - cos <p k ) ■ (2) 

fc=l 

The total electric energy T can be described as follows: 
T = \il c m( +\c g {v B -Vi ) 2 (3) 

fc=i ' ' 


where C g , V e and Vi denotes the capacitance of the gate 
capacitor, applied external voltage and the electric po¬ 
tential of node 1, respectively. Here, node 1 represents 
the superconducting island. 

Although the system Hamiltonian H has many energy 
levels, the system can be described as a two-level sys¬ 
tem (qubit) due to a strong anharmonicity by choosing 
suitable a. We show the a dependence of the energy of 
this system Fig. |l(b)| where E m (£ 12 ) denotes the en¬ 
ergy splitting between the ground (first excited) and the 
first excited (second excited) state. This clearly shows 
that system has an anharmonicity so that we can control 
only the ground state and first excited state by using 
frequency selectivity. 

| g) and \e) are the ground and the first excited state of 
the system Hamiltonian H = T + U for / = 0.5. In this 
regime, the ground state and the first excited state of 
this system contains a superposition of clockwise and an¬ 
ticlockwise persistent currents. Here, | L) = ^=(| g) + |e)) 
corresponds anticlockwise persistent current and |i?) = 
-^75(I e ) — Iff)) corresponds clockwise one. 

While / is around 0.5, due to the anharmonicity, we 
can consider only the ground state and first excited state 
in the Hamiltonian £, and so we can simplify the H into 
H ge spanned by | g) and \e) as follows: 

H ge = \ ( Acr z +e<J Y ) (4) 


where a z = |e)(e| - \g)(g\ and cry = -*|e)(g| + i\g)(e\ 
are Pauli matrices, A denotes the tunneling energy be¬ 
tween | L) and | R), e denotes the energy bias between 
| L) and | R). The energy of the qubit is described as 
£01 = Ve 2 + A 2 . 

In this paper, unless indicated otherwise, we fix param¬ 
eters as a = 0.2 and £,( 1 ) = 200 GHz and E j{k) /E c{k) 
ratio is 80. Here, £ c (fc) = e 2 /2Cj( k ) is charge energy 
of each Josephson junction. In this parameter regime, 
£01 is about three times larger than £12 as shown in 
Fig. 1(b) so that we could consider this system as an ef¬ 
fective two-level system. When / is set to be near 0.5, 
the derivative of the qubit energy against the magnetic 
flux l^fp-l takes the minimum value, so that the qubit 
should be well decoupled from flux noise, and we achieve 


¥>1 + V 2 + <P3 + <P4 = -27T/ 
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the maximum coherent times. We call this regime “opti¬ 
mal point”. On the other hand, we can control the value 
of £ by changing the value of /. When the energy bias e is 
much larger than the tunneling energy A, the persistent 
current states are the eigenvectors of the Hamiltonian so 
that we can read out the qubit state with SQUID[I] in 
{\L), 11?)} base. Here we show the dependence of s and 
A against magnetic field with no bias voltage applied in 
Fig. [2] It is worth mentioning that we can control the en- 



FIG. 2. The tunneling energy A and the energy bias £ against 
the magnetic flux f. e decreases monotonically as we increase 
/, while A is almost independent of /. 

ergy of the qubit by tuning the applied voltage V e while 
operating at the optimal point. We show the relationship 
between A and / with several values of V e in Fig. [3] 


III. ISING TYPE INTERACTION USING 
CAPACITIVE COUPLING 

A. Generating interaction between two-qubit 
system 

In this section, we show how to generate Ising type in¬ 
teraction using charge coupling for superconducting flux 
qubit. As a novel feature of our scheme, we use only 
external voltages to switch on and off the interaction be¬ 
tween two flux qubits. Unlike previous schemes, external 
magnetic field is not required to control the interaction 
in our scheme. Since the voltage can be applied locally 
compared with the case of applying magnetic field, we 
may have advantage in this scheme for scalability due to 
better individual addressability when we try to control 
individual qubits. 




(a) 


(b) 


FIG. 3. (a)The relationship between the external magnetic 
flux / and the energy of the qubit Eq\ with different voltage 
levels. Here, we set the gate capacitance Cg = 0.16 fF. (b) 
The relationship between / and Eoi with different a. Here, 
we set the gate capacitance Cg = 0.16 fF. 


Here, we show the circuit for our scheme using two su¬ 
perconducting flux qubits in Fig.|U The structure of each 


v (1 > v (2) 

e e 



FIG. 4. Two flux qubits 1, 2 are coupled via capacitance 
Cc^ 1,2 4 Each flux qubit is threaded by an external magnetic 
flux f( l \ and we can control the energy bias of the qubit via 
the magnetic flux. Node 1 and node 2 represent the super¬ 
conducting islands. JJ2 and JJ3 at each qubit have the same 
Josephson energies and capacitances that are a times larger 
than those of all remaining Josephson junctions. The electric 
potential of the island include node 1 (2) is ( Vf 2 ' 1 )- 


qubit is the same as that shown in Fig. |l(a)| When we 

apply an external voltage 14 on each qubit, the qubit 

(1 2 ) 

interact with each other across the capacitor Cc ■ We 
describe the details of this circuit in the following sub¬ 
sections. 
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B. Hamiltonian 


We now consider the electric energy and potential en¬ 
ergy of the circuit in Fig. [4] as follows: 

k= 1 1=1 x 7 

+lJ2 c 9 l) (vP-vV ) 2 

1=1 


+\cp»(vp >-if>) 2 (5) 

< 6 > 

fc=l 1=1 

2 

= T + U = E h a + H b (7) 

1=1 

h a = \Y.cf m fel ) ) 2 +\cf (yP - vp) 2 

k=l ' 7 

+ Z) s ?(fc) ( 1-cos ^°) ( g ) 

k=1 

Hb = \cm (v; (1) -f/ 2) ) 2 (9) 


where Cg Z \ 14^', and denotes gate capacitance, 
external magnetic flux, applied external voltage, and the 
electric potential of the island including node l for the l 
th qubit respectively. Here, node l represents the super¬ 
conducting islands. 

For an arbitrary /, we can derive the effective four- 
level Hamiltonian H ge of the eigenspace spanned by \gi) 
and | ei) from H to t a i■ Here, | gi) and | ej) correspond to 
the ground state and first excited state of the I th qubit 
without interactions for /W = /^ == 0.5. We expand 
Htotai by \g)i and \e)i. The effective Hamiltonian H ge 
becomes as follows: 

4e= E \viV 2 )(VlV2\H to t a i\viV2)(viV 2 \ 

vi€(gi ,ei) , 1 ) 26(92 ,e 2 ) 

= E ( a ( iv * +£<,) <) + s ^ 2 ' (io > 

1 = 1 z 


where g denotes the Ising type interaction strength be¬ 
tween qubit 1 and 2. We show the change of the qubit 
energy Eq\ and the interaction strength g as a function of 
applied voltages in Fig. [5] Large interaction strength and 
small derivative of qubit energy against voltage can be 
achieved by the large coupling capacitance Cc between 
each qubits. This seems to show that one can suppress 
errors by increasing Cc. We discuss about the errors dur¬ 
ing controlled-phase gate operation in following section. 



FIG. 5. The voltage dependence of the qubit energy A and 
the interaction strength g between two qubits of the circuit 
in Fig. [4] Here, both of the gate capacitance C l g = 0.077 fF, 
coupling capacitance C c = 0.077 fF. 


C. Effects on interaction from change in electric 
field 


To evaluate the performance of our scheme, we focus 
on two types of errors. Firstly, we analyze the dephasing 
errors due to the fluctuations of applied voltage. We 
define this type of error Cd and dephasing time T 2 as 
follows: 


trr) 7T” 

= 77V, t cp = —, X 2 = 

T 2 4 g 


^ 21 1 Sv 


( 11 ) 


dv 


where we assume t cp « T 2 . Here, t cp denotes the neces¬ 
sary time to perform a controlled-phase gate with Ising 
type interaction, v denotes the external voltage of each 
qubit, and Sv denotes the fluctuation width of v. It is 
worth mentioning that td has a linear relationship with 
Sv. To make 6d smaller, We should obtain a parame¬ 
ter set where the absolute value of the gradient of the 
qubit energy Eq 1 is small and the interaction strength g 
is large. 

Secondly, we investigate the jitter error of a two-qubit 
gate operation. The Ising type interaction can implement 
the controlled-phase gate 


Uc£ ] (t) = exp 




( 12 ) 


where g denotes the interaction strength in Eq. m, 
t = 4L denotes the time to apply voltages, and Uq’ z ' 
denotes a controlled-phase gate between qubit 1 and 2. 
By performing the controlled-phase gate on two qubits 
which are initialized to |++)i 2 state, we can obtain the 
two-qubit cluster state. But, the applied voltages may 
not create the desired state due to error in the timing 
t' = t + St , where St is timing jitter. We introduce the 
controlled-phase gate U^ z ' > (f) including the timing error 
to calculate a gate fidelity Fez = |(</ ) | < / ), )| 2 with 

10) = Ucz\t)\++), \4>) = Ucz\t') |++). (13) 
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Here, we define the timing error e t , m = 1 — Fez, and the 
local error q oc (= + e t j m ). We show the e; oc against 
the applied voltage V e with the particular values of Cc 
in Fig. [G] The threshold of local errors for fault-tolerant 



FIG. 6. The total local error ei oc {= £d + ctim) as a function of 
voltage with different coupling capacitance Cc. Here, we set 
the fluctuation width of voltage Sv = 0.21 gV and the timing 
jitter St = 50 psec. Dashed line denotes an error of 0.1%. 

quantum computation is known to be around 1%. Also, 
it is known that, if the error rate is close to the threshold, 
the necessary number of qubits for the computation dras¬ 
tically increases [4^, |4l|. Therefore, we set the threshold 
to ei oc = 0.1%. As shown in Fig. 0 we can increase the 
coupling strength g by increasing Cc. Meanwhile, the 
strong coupling strength causes the large timing error. 
Therefore, as shown in Fig. [G] the optimal voltage exists 
for each of the Cc which minimizes the total local error. 
In addition, by increasing Cc, the total error tends to be 
smaller. This result shows that the large Cc has an ad¬ 
vantage for quantum error correction against local errors. 
However, for nmlti-qubit systems, increasing Cc causes 
a different problem. Unwanted interaction strength be¬ 
tween non-nearest neighbor qubits increases due to the 
large Cc. For this reason, the Cc should be set to be 
around 0.075 fF. The detail of this will be discussed in 
Section IIVBI 


IV. MULTI-QUBIT SYSTEM 

In this section, we generalize our scheme to multi-qubit 
system. Firstly, we discuss how to control the capaci¬ 
tive interactions between superconducting flux qubits via 
applied voltage. Secondly, we show how to apply our 


scheme to generate a two dimensional cluster state using 
superconducting flux qubits arranged on square lattice. 

A. Generating interaction between multi-qubits 

system 

Here, we discuss the interactions between capacitively 
coupled N flux qubits that are arranged in one dimen¬ 
sional line as shown in Fig. [7] For simplicity, we assume 
homogeneous flux qubits. denotes the external mag- 


y(D y(2) y(3) y(N) 

e e e e 



FIG. 7. A flux qubit at the site j(l < j < N) couples with 
the nearest neighbor qubits via capacitance C'c ,:J±1 ' ) . For sim¬ 
plicity, we assume homogeneous flux qubits. Each node j 
represents the superconducting islands. Each qubit has four 
Josephson junctions. Two Josephson junctions directly con¬ 
nected to the node have the Josephson energies and capaci¬ 
tances that are a times larger than the other two Josephson 
junctions. 

netic flux through the loop of the j th qubit. When all 
flux fti) are 0.5, the system Hamiltonian is described as 
follows. 

N N 

H = £ 2 A (i) a{ z + J2 9(\i-i'\)Vz )<J z ] ( 14 ) 

z=i z i,i'=1 

where Am denotes the energy of the l th qubit, 
denotes the interaction strength between each pair of 
qubits at a site [1,1'), and \l — l'\ denotes the site dis¬ 
tance between these qubits (e.g. when qubit l and l' are 
nearest neighbor pair, |/ — l'\ = 1 .). 

B. Generation of a one dimensional cluster state 

Non-nearest neighbor interactions cause spatially- 
correlated errors that are difficult to correct by quan¬ 
tum error correction. In this subsection, we show the 
way to evaluate this error. We define the ratio between 
nearest neighbor interaction g[= <?(i)) and next-nearest 

neighbor interaction g^ 2 ) as R (= w h ere all qubits 

are applied voltage V e . We show that the interaction 
strength g[\l — l'\) decreases exponentially as the site dis¬ 
tance — V | increases, and the Ratio R depends on the 
coupling capacitance Cc between each qubit. This is a 
striking feature in our scheme using voltage for the con¬ 
trol of the qubit-interaction, because the effect from any 
control lines can decrease only polynomially against the 
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site distance if one uses magnetic field for the control. 
We show the interaction strengths of 6 qubits system as 
a function of Cc in Fig. [8] 



FIG. 8. The Cc dependence of the interaction strengths and 
the coupling ratio R where g(\l — l '\) denotes the 

interaction strength between a pair of qubits at a site (l — l 1 ). 

If we apply voltage on all qubits, interaction occurs 
between such qubits. The total error enon caused by 
non-nearest neighbor interactions on j th qubit during 
controlled-phase operation is calculated as follows: 

N/2 N/2 

4?on = 9(n)tc P m ( n) = ^ ^-R ( " _1) m (n ) (15) 

n =2 n =2 

where n denotes the site distance between the j th qubit 
and the coupled non-nearest neighbor qubits, de¬ 

notes the number of such non-nearest qubits. 

Such the existence of the spartially-correlated error will 
increase the threshold for quantum error correction [42} . 
Large capacitance tends to decrease local errors as shown 
in Fig.[6l while large capacitance induces more spartially- 
correlated errors as shown in Fig. [8] However, when we 
consider the spatially-correlated error, the error thresh¬ 
old value of the surface code is not well studied. Thus, 
we set the upper bound of the spatially-correlated error 
on each qubit c non < 10 q 00 which is an order of mag¬ 
nitude smaller than the threshold of local error for sur¬ 
face coding scheme. If this condition is satisfied, we as¬ 
sume that spatially-correlated error is enough to perform 
a fault-tolerant quantum computation. When we apply 
voltage on all qubits to perform controlled-phase gates to 
all pairs of nearest neighbor qubit, a range of values that 
the coupling capacitance Cc can take while satisfying the 


above condition is smaller than the proper range of Cc 
discussed in Subsection IIII Cl Therefore, we do not apply 
voltage on all qubits but apply voltage on some of them. 
We choose pairs of nearest neighbor qubits that we will 
apply the voltage, and we set a site distance p between 
the pairs. Then, if R is small enough, e non of each qubit 
is the following equation: 


, 0 ) 


N/2 


y -r^ 


n—p 


1 


- 10000 


(16) 


where p is the site distance between qubits applied by 
voltage. 

Since there are many parameters on the interaction 
Hamiltonian, it is difficult to find an optimum set of 
parameters that minimize both of local and spatially- 
correlated errors. Therefore, we fix the following param¬ 
eters: a = 0.2, S v = 0.21 pV, S t = 50 psec. To determine 
a minimum site distance p while suppressing the corre¬ 
lated errors to be under 0.01 %, we show the Cc and V 
dependences of the errors with p = 4, 5 in Fig. [9] As 



FIG. 9. (a)The Cc dependence of the correlated errors. 
Dashed line corresponds to an error of 0.01 %. (b)The V 
dependence of the total errors. Dashed line corresponds to an 
error of 0.1 %. 


shown in Fig. |9(a)| when p = 4, the e non exceeds 0.01 
% around Cc = 0.04 fF. We cannot sufficiently suppress 
local errors using coupling capacitance smaller than 0.07 
as shown in Fig. [6] Thus, the site distance p should be 
larger than 5. Meanwhile, when p = 5, the e non exceeds 
0.01 % around Cc = 0.09 fF. Then the total error of the 
controlled-phase operation can be sufficiently suppressed 
to be less than 0.1 % using the coupling capacitance Cc 
around 0.077 fF as shown in Fig. |9(b)| Therefore, it is 
preferable that the site distance p = 5 be selected. In or¬ 
der to adopt sufficiently large coupling capacitance such 
that the ei oc below 0.1 %, we need to choose sufficiently 
large p such that the e non below 0.01 %. We discuss about 
the way which can further reduce p in the following. 

The p determines the maximum number of controlled- 
phase gates that are performed simultaneously on the 
same system. For example, we can perform L^yrJ + 1 
controlled-phase gates in parallel using iV-qubits one di¬ 
mensional system. If we can use the smaller p without 
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adding extra errors, we can perform more controlled- 
phase gates in parallel. So that we can generate a cluster 
state within a shorter operating time. For this purpose, 
we introduce the spin echo technique where implementa¬ 
tion of a 7r pulse (single qubit ax rotation) to the target 
qubit could refocus the dynamics of the spin so that ef¬ 
fects of interactions on the target qubit should be can¬ 
celled out. We apply two tt pulses to pairs of qubits to 
suppress spatially-correlated errors. For example, we set 
three qubits in a raw and apply voltage vj to the n 
th qubits {n = 1,2,3) as shown in Fig. [TOl where Ve 

(2) (3) 

and Ve are equal, V e is an arbitrary voltage, and the 
strength of interaction between qubit 1 and 2 is g. We 


e 2 e ' 3 
6(1) -.6 (1) 


Q: O :0 

S'(2) 

k pulse t 

.O 

<> t 

o—o o 


kH + + +>i; 


1 - 2 5 12 “ *cd 2 


I. A -iH^t -iH u t 

\y/)-e 12 e 23 e 13 


2 

+ + A 


\V') = e 
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+ + + 


( 1 , 2 ) I 


FIG. 10. When we perform a n pulse on qubit 1 and 2 at t = 
t C p/2 , the nearest neighbor interaction between qubit 2 and 3 
and the non-nearest neighbor interaction between qubit 1 and 
3 are cancelled out. In such way, we can perform controlled- 
phase gate without changing the state of other qubits. 


set each qubit to be prepared in |+) state, let the state 
evolve for a time t cp /2 , perform two tt pulses to qubit 1 
and 2, and let the state evolve for a time t cp / 2. The final 
state become as follows: 


U\+ + T) 123 = (jTO) 12 T |— 1 ) 12 ) ® |T)3- (17) 


Here, the interactions and 3(2)0’k (T z^ are 

cancelled out due to the tt pulses and we obtain a cluster 
state between qubit 1 and 2. 

This method can be applied with the case of arbitrary 
number of qubits. The general rules are follows: let us 
consider a pair of qubits. If we perform tt pulses on both 
of qubits, the interaction between them is not affected by 
these pulses. On the other hand, if we perform tt pulse 
on one of them, the interaction between them is cancelled 
out. These properties would be crucial for generating a 
cluster state as we will describe. 

For generating a large one dimensional cluster state 
using N qubits of the circuit in Fig. [3 we show the pro¬ 
cedure as follows: 


Step 1: We apply voltage to (3n — 2) th and (3n — 1) 
th qubit for performing controllcd-phase gates be¬ 
tween (3n — 2) th and (3n — 1) th qubit where 


n = l,2,..,L^J. 


vv vv vv vv 

v e(l) v e(2) v e(4) v e<5) v e(7) v „(8) v e(10) v e(ll) 



FIG. 11. The 3-step procedure for generating a one dimen¬ 
sional cluster state. Step 1. We initialize 3n — 2 th and 3n— 1 
th qubits in |+). Here, n = 1,2, ••, [ Af ^ 1 J where [zj is the in¬ 
teger part of x. After that we apply voltage on 3n — 2 th and 
3n — 1 th qubits. Let the state evolve for a time perform 
tt pulses to 6n — 2 th and 6n — 1 th qubits, and let the state 
evolve for a time ^. After these operations, controlled-phase 
gates have been performed between qubit 3n — 2 and 3n — 1. 
Step 2. We initialize 3n th qubits in |+). After that, simi¬ 
lar to the Step 1, we perform controlled-phase gates between 
qubit 3n — 1 and 3n. Step 3. We initialize 3n + 1 th qubits 
in |+). After that, similar to the Step 1 and 2, we perform 
controllcd-phase gates between qubit 3n and 3n + 1. 


Step 2 : We apply voltage to (3n — 1) th and 3n th 
qubit for performing controlled-phase gates be¬ 
tween (3 n — 1) th and 3n th qubit where n = 


1 , 2 ,", AAJ- 


Step 3: We apply voltage to 3n th and (3n + 1) th 
qubit for performing controlled-phase gates be¬ 
tween (3n — 1) th and 3n th qubit where n = 

1,2,", L^J- 

At each step of the above procedure, LkkJ controlled- 
phase gate are performed in parallel. At each step, 
we will perform the following procedure to perform the 
controlled-phase gate. Firstly, prepare the qubit state in 
| T) • Secondly, let the state evolve for a time t = ^ ac¬ 
cording to the Hamiltonian described in Eq. Ildl Thirdly, 
perform the n pulses to suppress the non-local interac¬ 
tion. Finally, let the state evolve for a time t = t cp - 
We show the details of these operations in Fig. [TT] and 
explain how the non-local interaction is suppressed in 
Fig. Q2] When all coupling capacitance are Cc < 0.077 
fF, the spatially-correlated error on each qubits become 
as follows: 


N/2 


d=E^ (n " 1)m 

n =5 


r n ) ~ —(R 4 + 2R 5 ) < - 

(n) 4 y 10000 


.(18) 


The k th qubit is affected by mainly three non-local 
interactions as shown in Fig. [12] The strength of the 
largest interaction is gR 4 , and the strength of the other 
two interactions are gR 5 . The remaining non-local inter¬ 
actions are negligibly small. 










Voltage Voltage + n pulse 


a T?5 (£-6) (*) 

g ( l)K <7 Z (7 


Voltase Voltage + n pulse 
_*/ V 



gn)R °z 


(k)(k+6) 


FIG. 12. The influence of non-local interactions. During the 
controlled-phase gate, each target qubit are affected by non¬ 
local interactions. We show the strength of mainly three non¬ 
local interactions with k th qubit. These interactions are not 
cancelled out by n pulses. 


C. Generation of a two dimensional cluster state 

Next, we show how to generate a two dimensional clus¬ 
ter state using TV 2 flux qubits arranged on TV x TV square 
lattice. We show a part of the circuit in Fig. [131 fV' k ) 


y( m , 


1 , 111 + 1 ) 


y(m 


+l,m+l) 



FIG. 13. Physical circuit for generating a two dimensional 
cluster state. These four qubits correspond to the qubits 
surrounded by dot line in Fig. 1141 Two Josephson junc¬ 
tions directly connected to a node (the superconducting is¬ 
lands) have the Josephson energies and capacitances that are 
a times larger than the other two Josephson junctions. Every 
flux qubit at site (j,k) couples with the four nearest neighbor 
qubits via capacitance C'c^’ fe ^ ±1 ' fc±1 ^. 


denotes the external magnetic flux through the loop of 
the qubit at site (j. k). Here, (j. k) corresponds to the 
lattice point. When all flux fU’ k > are 0.5, the system 
Hamiltonian is described as follows: 


k= E 

( l,m ) 


(i , m) 

2 z 


V'' (Z,m) (l' ,m') / 1n \ 

+ 2 ^ 9(\l-l'\+\m-m'\)^Z °Z ( 19 ) 

((Z,m),(Z / ,ra / )) 


where A denotes the energy of the qubit at site 
(Z,m), gnt-i'\+\m-m.'\) denotes the interaction strength 


between each pair of qubits at site ( l,m ) and (Z',m'), 
and \l — l'\ + \m — m'\ denotes the site distance between 
these qubits. 



FIG. 14. Schematic of our procedure for generating a two di¬ 
mensional cluster state by graph state representation. Circles 
correspond to qubits, dashed lines correspond to electrically 
connection via a capacitance, solid-lines correspond to en¬ 
tanglement between qubits, and numbers show the order in 
which controlled-phase gates are performed by our procedure. 
White circles denote separable qubit, and gray circles denote 
qubits constituent of cluster state(s). 

Here, we show the 12-step procedure as follows for gen¬ 
erating a two dimensional cluster state. 

Step 1-3: We perform (TV— l)|_xJ controlled-phase gate 
to generate [-yj one dimensional cluster states us¬ 
ing qubits located in the 4m—3 (to = 1, 2, L^y^J) 
row in the same way as shown in Fig. 1111 Then the 
spatially-correlated error of each qubit in the 4m—3 
row is smaller than ln ^ nn . We show the outline of 
these steps in Fig. fl4fa). 

Step 4-6: We perform (TV — l)L^y-^J controlled-phase 
gate to generate one dimensional cluster 

states using qubits located in the 4p — \{p = 
1, 2, (A+iJ) row in the same way as above. We 
show the outline of these steps in Fig. Qjjb). 

Step 7-9: We perform (TV— l)[ ; yJ controlled-phase gate 
to generate a two dimensional graph state as shown 
in Fig. Etc) using qubits located in the 4m — 3 
column across [-yj one dimensional cluster states. 
We show the outline of these steps in Fig. fl4T c). 

Step 10-12: We perform (TV— 1) L^-y^J controlled-phase 
gate to generate a two dimensional cluster states 
using qubits located in the 4p — 1 column. We show 
the outline of these steps in Fig, fill' d). 
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FIG. 15. Operations and the influence of non-local interac¬ 
tions in generating a two dimensional cluster state. In this 
step, we apply voltage to qubit at site (3n — 2, 4m — 3) and 
(3n — 1,4m — 3). Let the state evolve for a time -tt, perform 
7r pulses to qubit at site (6n' — 5, 8 m! — 7), (6n ' — 4,8 m! — 7), 
(6n'— 2, 8m 1 — 3), and (6n / —1, 8m' —3), and let the state evolve 
for a time So that controlled-phase gates can be imple¬ 
mented between the pair of qubits at site (3n — 2, 4m — 3) and 
(3n—1, 4m—3). Here, m = 1,2, », [^J, m' = 1,2,-, 
n = 1,2, ••, [ ^j" 1 J, and n' = 1, 2, ••, ]. Each target qubit 

is affected by non-local interactions from qubits on the same 
row and other rows. We show mainly five non-local interac¬ 
tions with the qubit at site (j,k). These interactions are not 
cancelled out by 7r pulse. 


We show the details of each step of above procedure 
for generating a two dimensional cluster state in Fig. 1151 


During each step, a part of the non-local interactions 
are not cancelled out by 7r pulses. When all coupling 
capacitance are Cc < 0.077 fF, the spatially-correlated 
error on each qubits become as follows: 

N/2 

= .( 20 ) 

71=5 

The qubit at site (j, k) is affected by mainly five non¬ 
local interactions as shown in Fig. [15] The strength of 
the largest interaction is g^R 1 2 3 4 , and the strength of the 
other four interactions are g^R 5 6 7 8 9 10 . The remaining non¬ 
local interactions are negligibly small. 

V. CONCLUSION 

In conclusion, we suggest a new way to generate Ising 
interaction between capacitively-coupled superconduct¬ 
ing flux qubits by using an applied voltage, and we also 
show architecture about how to make a two-dimensional 
cluster state in this coupling scheme. Unlike the stan¬ 
dard schemes, our scheme does not require to change the 
applied magnetic field on the flux qubit for the control 
of the interaction. Since applying local voltages is typi¬ 
cally much easier than applying local magnetic flux, the 
scheme described in this paper may have advantage to 
perform two-qubit gates on target qubits without affect¬ 
ing any other qubits. Our result paves the way for scal¬ 
able quantum computation with superconducting flux 
qubits. 
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